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Abstract 

The measurement of the forward- backward asymmetry at LHC could 
be an important instrument to pinpoint the features of extra neutral gauge 
particles obtained by an extension of the gauge symmetry group of the 
standard model. For definitiveness, in this work we consider an extension 
of the gauge group of the Minimal Supersymmetric Standard Model by 
an extra anomalous U(l) gauge symmetry (MiAUMSSM). We focus on 
pp —¥ e + e~ at LHC and use four different definitions of the asymmetry 
obtained implementing four different cuts on the directions and momenta 
of the final states of our process of interest. The calculations are performed 
without imposing constraints on the charges of the extra Z's of our model, 
since the anomaly is cancelled by a Green-Schwarz type mechanism. Our 
final result is a fit of our data with a polynomial in the charges from which 
to extract the values of the charges given the experimental result. 
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1 Introduction 



One of the most motivated extensions, from a theoretical point of view, of the 
Standard Model (SM) and Minimal Supersymmetric Standard Model (MSSM) 
of particle physics is obtained by enlarging the gauge group of the theory by 
admitting extra [/(l)'s. Such extensions are natural at low energy for models 
coming from grand unified theories and string theories (see [T] for a recent 
review). In the string inspired scenarios the anomalies of the extra £7(1) 's are 
cancelled by the Green-Schwarz (GS) mechanism. To explore such possibility 
we will use an extension of the MSSM which from now on will be dubbed 
MiAUMSSM. An alternative version of this model which admits spontaneous 
supersymmetry breaking was also formulated in [2] , but in this work we will use 
the original formulation of [3] . The phenomenology of the MiAUMSSM has been 
investigated in different directions. Assuming that the lighest supersymmetric 
particle (LSP), a candidate for dark matter, comes from the anomalous sector 
of the model [H [5] , the relic density of such LSP was computed and proved to 
be compatible with the experimental data of WMAP [5] . Furthermore in [7] the 
decays of the next to lightest supersymmetric particle (NLSP) into the LSP has 
been considered while in [B] the features of a possible signature of the model at 
LHC has been considered by concentrating on a particular particular radiative 
decay of the NLSP. 

In this paper we will further develop the phenomenology of the MiAUMSSM 
by computing the forward-backward asymmetry which is induced in the final 
states of the process pp — > e + e~ by keeping into account the new gauge boson, 
Z' , associated to the extra [/(l)gauge symmetry. The couplings (charges) of 
this particle to the others present in our model are not fixed by the requirement 
of gauge anomaly cancellation and can be determined only by experiment. Our 
aim is to show that such measurement is feasible and that it can distinguish 
among the different possible scenarios [S]. Since at LHC the colliding beams are 
made of the same particle, to generate an asymmetry in the final state, some 
cuts on the parameter space have to be necessarily performed. Each possible 
cut leads to a different definition of the asymmetry. In this work we will use 
four different sets of cuts to show that our results are not dependent from these 
choices. 

This work is organized as follow: in sec. [2] we briefly review the main 
features of the model which we are going to use. In sec. [3] we will discuss 
the four different definitions of the asymmetry we will use, in sec. [4] we will 
describe our calculations and collect the results which are finally discussed in 
the conclusions. 

2 Model definition 

Our model [3] is an extension of the MSSM with an extra U(l). The charges of 
the matter fields with respect to the symmetry groups are given in table 1. 
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Table 1: Charge assignment. 



The gauge invariance of the model implies: 

Qu" = —Qq - Qh u 
Qd<= = —Qq + Qh u 

Qe< = ~Ql+Qh u (1) 
Qn d = —Qh u 

So there are only three free charges introduced by the extra symmetry: we can 
choose Qq, Ql and Qh u without loosing generality. The anomalies induced 
by this extension are cancelled by the GS mechanism: there are no further 
constraints on the charges. 

To evaluate the asymmetry associated to the full process pp — > e + e _ we have 
performed the calculation of the cross section of the subprocess qq — > e + e~, 
which we report in Appendix[Aj In the next section we evaluate the convolution 
of this differential cross section for the specific definitions of asymmetry we will 
adopt. We take the mass of our Z' to be 1.5 TeV. There are two main reasons 
for this choice: on the one hand we wanted a sizeble Z' production (see [3], 
where there are results for a Z' mass of 1 TeV). On the other hand this mass 
value allows a comparison with the results contained in literature |10j . Anyhow 
our analysis could be repeated for arbitrary value of the Z' mass. 



3 Asymmetry definition 

Due to the fact that the initial pp state is symmetric, the asymmetry at LHC 
is zero if we integrate over the whole parameter space. However the partonic 
subprocess qq — > e + e~ is asymmetric. We can keep this asymmetry by imposing 
kinematical cuts, which are anyway inevitable because of the limits imposed by 
the detector. There are many possibilities to perform these cuts and each of 
them leads to a different definition of the asymmetry. In this work we have used 
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the four definitions of the asymmetry which are collected in [10] : 
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where a is the total cross section after integrating with the partonic distribution 
functions (PDFs). 

The first two asymmetries are defined in the center of mass (CM) frame. 
The forward-backward asymmetry Are r [JJ El [UJ H21 1131 E] has a cut on the 
rapidity Y of the e~/e + pair, while the one side asymmetry Ao |15[ 1 1 6 j has 
a cut on p ze - e +, the total momentum associated to the final states (e~e + ) 
moving longitudinally along the beam direction chosen to be the z axis. In the 
following, for the sake of simplicity, we will drop the subfix in the definition of 
the rapidity 
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In ( 10 ) this rapidity will be expressed in the CM in terms of the partonic vari- 
ables Xi t 2- E e - e + is the sum of the energies associated to the two particles. The 
other two asymmetries are defined in the laboratory (Lab) frame. The variable 
Y e ± is the pseudo-rapidity associated to the single particle e ± and expressed as 



Y„. 



log (^tan(6» e 



72) 
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with 6 e the angle of the outgoing fermion with respect to the z axis. In 
this case the kinematical cut is over the rapidity in the Lab frame which is 
denoted by Yq and which will be introduced later on. The central asymmetry 
Ac [III HI HS1 EQl [21] is calculated integrating in the angular region centered on 
the axis orthogonal to the beam, while the edge asymmetry Ae [22] is defined 
in the complementary region. 

The explicit expression of eq. Q in the CM frame is: 
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the forward and backward contributions, respectively. f q /q{xi) are the PDFs of 
q/q. C cu t is the domain of integration, that depends on the type of asymmetry 
that we want to calculate according to the definitions (2|5). 
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The couples of variables (an, x 2 ) and (s, Y) are not independent. In fact, their 
definition is: 

a = Sana* Y=hog(^-) (10) 
2 \x 2 / 

where S = (14 TeU) 2 is the total squared energy of the accelerator. 

Since we have calculated the cross section of the process (see Appendix |A|) 
that we are going to study with respect to s, we perform the change of variables 
(an, x 2 ) — > (s, Y) in the expression pi). The Jacobian of this transforma- 
tion is J = 1/S. Now we focus on the integral j_ Ycut dY: if we perform the 



change of variable Y — > — Y", because of (10), this corresponds to the exchange 
an <-> a* and consequently to the exchange of forward with backward (F <-> B). 
Summarizing 

/•«> 

dy^/ g (a; 1 )/ g -(ar 2 )(F± J B)= / dr^/ 9 (x 2 )/ (7 (x 1 )(±F+ J B) (11) 
-00 g Jy c „ t q 

Using this result we obtain the formula that we will use in this paper, that is: 

JdsJj+2dYj: q f- s (x 1 (s,Y),x 2 ( S ,Y))(F-B) 

AftFB = r—, (12) 

/ ds J /+~ dY J2 q Y), x 2 (s, Y))(F + B) 

with 

fqq{x 1 ,X 2 ) = fq(x 1 )fq{x 2 ) ± fq(x 1 )f q (x 2 ) (13) 

The formula for is very similar and it is obtained by replacing Y e - e + with 
p z . This leads to a different cut and a different Jacobian. 

The expressions for Ac and A# are also very similar between them and 
we present them together. In these two cases the cut is not performed on the 
rapidity or on the longitudinal momentum but on the angle 9 cut (therefore on 
the limits of integration for F and B). These asymmetries are defined in the Lab 
frame, because the Lorentz transformation from the CM frame "squezees "the 
final particles [17] • The angles of the outgoing e T with respect to the z axis, 
denoted by ±9 respectively in the CM frame, are replaced by 9 e and 9 e in the 
Lab frame where the outgoing e ± no longer have the same direction. Therefore 
in the Lab frame, instead of the definitions (j9J), we have: 

(F/B) c = *x dcoB ge ( - /+) (U) 

L.UL dcOSP e 

(F/B) E = f- CUt da_ dwa0 e^+-> + f + 1 ^_dcOS0 e< - /+) 

V ' /fc J -! dcos0 e(_/+) ds Jcut dcos0 eC /+ 'ds 

where the extremum of integration cut is defined in terms of 9 cut as cut — 
cos 6< cut . Then Q and Q become: 

A c /e = J ds J^ s dx 1 J^2f q (x 1 )f q (x 2 = 1^)( f c/e- B c/E ) (15) 

J is now the Jacobian of the transformation (an, a*) ~~ ^ (xi,s). 
The cut parameter which enters in the definitions Q and ^ is not directly 
Qcut assoc i a ted pseudorapidity = — log (^tan(6> CMt /2)^ . For further 

details on the calculations sketched in this section, see [25] , 
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4 Asymmetry calculation 



In this work we have calculated the asymmetry in two different ways. First 
we have used a numerical code that we have written using Mathematica. This 
code uses the cross-section computed in Appendix [A] to numerically compute 
the integrals discussed in the previous section. As a second check we have 
repeated the same computation using the HERWIG package [HI [55] , that we 
have modified to calculate the asymmetry. We have chosen to repeat twice our 
computation for two main reasons: the first one is that in this way we can have a 
cross check between our results; the second is that these methods have different 
peculiarities that we want to use. For example, the numerical integration is less 
computer time consuming for the Mathematica code, which helps in establishing 
the dependence of the asymmetry from the free charges of the model. At the 
same time the HERWIG package permits to study how the cuts influence the 
rate of production of our final state. For these reasons we have performed the 
basic calculation (i.e: the asymmetry optimization) using both methods. We 
remark that all the results that we will show are strongly dependent on the set 
of PDFs used to calculate them and that this leads to a systematical error. In 
the following we do not show results for different sets of PDFs. For what the 
statistical error is concerned, we have estimated it using the formula |10j : 



err = \/ — — — ~ . — (16) 

where N F / B are the forward/backward events, N is the total number of events 
and C is the luminosity. In the following we show the estimated errors for the 
asymmetry definitions keeping L = 100 

We aim to use the asymmetry to distinguish our model from the MSSM or other 
models which include an extra U(l). As an example, in fig [I] we plot Arfb for 
our model and for the MSSM. It is evident what is the effect of our Z' resonance 
on the asymmetry, that is the presence of a dip in the trend of the asymmetry 
around the value of mass of the Z' . If a model has a different Z' resonance 
we also expect differences in the shape of the curve showed in fig. [I] In the 

0.8 r 




Figure 1: FB Asymmetry in MiAUMSSM and MSSM versus M e+e - around the 
resonance of the Z' . y cut is set to 0.4, while the anomalous charges are fixed to 
the sample values Ql = 1, Qh u =1/2 and Qq = 3/4 

following we will perform the asymmetry calculation around the peak region, 
that is for Mz' — 3Fz> < M e + e - < Mz> + 3Tz>, where Yz> is the total decay rate 
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of the Z' . As we previously remarked, this determines the integration domain, 
that is {Mz> - 3IV) 2 < s < [M z < + 3IV) 2 . We also compare our results with 
the ones obtained for the Sequential Standard Model (SSM), in which there is 
an extra Z' boson which has the same couplings to fermions such as as the SM 
Z boson [3[2S]. In fact, in this model, the extra U(l) gauge symmetry is defined 
by the Lagrangian: 

gz. J%,Z' M = -J- $ff (gl - oils) TpfZ'^ (17) 
lc w f 

where c w is the cosine of the Weinberg angle and g, g[,g[ are the vector and 
axial coupling costant of the Zq boson of the SM. See table [3] for further details 
on the values of the quantum numbers. 

4.1 Optimized asymmetry 

As shown in [10] the asymmetry magnitude is not a good function to optimize. 
A better choice is, instead, the statistical significance: 

Sig = AVC^ (18) 

where A can be any of the previously defined asymmetries, C is the LHC inte- 
grated luminosity, which we take to be 100 fb^ 1 . 

In figure [2j [3] and [4] we show the results for the on-peak significance of the Mi- 
AUMSSM and SSM for all the definitions of asymmetry that we use. There is 
a good agreement between the results obtained by using the Mathematica code 
(Figs[2]and[3| and those obtained using the event simulator HERWIG (Figj4j . 
So we are confident that our results are reliable also when we will use them to 
calculate other observables, e.g. the dependence from the charges of the asym- 
metries and the significancies. The best cuts are those that maximize the 
significance. For the SSM, we find the same values as in [10]. We list the best 
cuts of the MiAUMSSM in table [2] As in [10] . we expect that the best cuts are 
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Table 2: Best cuts for the on-peak e + e asymmetries 



nearly independent from the charges and depend only on the Z' mass and the 
partonic distribution functions. Moreover they are also essentially independent 
from the specific model chosen as it is confirmed by our analysis. As a further 
check we have performed simulations with the SSM. We used the same settings 
of [10] . obtaining very similar results for all the cuts, confirming the reliability 
of our numerical codes. We used the SSM not only for having a check of the 
validity of our results, but also to have results that can be compared with those 
of the MiAUMSSM. 

4.2 Dependence on the charges 

Now we want to use the best cuts previously found to study the asymmetry 
in function of the free charges of our model. We have studied the value of 
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Figure 2: Significance as a function of the corresponding cut for the four def- 
initions of asymmetry for on-peak events in the MiAUMSSM (computed with 
Mathematica) 



the four asymmetries keeping alternatively one of the charges fixed to and 
varying the others two from —1 to 1. We choose these ranges because in the 
SM all the charges are of this order. Furthermore in [5] we have found that 
— 1 < Qh u ;$ 1 for the model to be consistent with the WMAP data on dark 
matter. Then, out of simplicity we have used the same region for the other two 
charges. As said, we have obtained the following plots using Mathematica to 
perform the numerical integration. Some of the results we have obtained are 
showed from figure [5] to figure [lOj In these plots the darkest color areas are 
those with the lowest absolute values of the asymmetry while the greatest values 
lie in the lightest color region. In addition, these plots are almost symmetric 
under exchanges in the signs of the charges. The contour plots with Qh u = 
(figures [9] and 10) are almost invariant under the exchange Qt — > —Qi of the two 
remaining charges. Those with Ql = or Qq = are almost symmetric only 
under the change of sign of both the two unfixed charges. So the asymmetry as 
a function of the charges must reflect these sort of symmetries in its polynomial 
dependence on the charges. These implies that if we try to fit the asymmetry 
with a rational function (that is the best choice, given the definitions |2][5])) we 
will have constraints on the coefficients of the fit. 



4.3 Number of events 

We already mentioned that to obtain a non zero asymmetry at LHC we have 
to impose cuts in the parameter space. Obviously these cuts will diminish the 
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Figure 3: Significance as a function of the corresponding cut for the four defini- 
tions of asymmetry for on-peak events in the SSM (computed with Mathemat- 
ica) 



number of events that we can use to measure the asymmetry, so it is important 
to be sure that they do not drastically affect the set of data we have at our 
disposal. To study the ratio between the number of events obtained applying 
the cuts and the total number of events expected in our channel of interest 
(pp — > e + e~) we have utilized the HERWIG package. We have studied the 
ratio Ni /N to ti , where Ni is the sum of the forward and backward events for 
the i-th definition of asymmetry and Ntoti 1S ^ ne number of events that we have 
generated with HERWIG. We have performed the calculation of Ni/N to ti in two 
cases: 

• on peak invariant mass, variable cuts 

• variable invariant mass, fixed cuts 



The results are showed in figures 11 and 12 We can see that after the im- 
plementation of the cuts we are left with the 65 — 75% of the total number of 
events for the FB and O asymmetry, while for the C and E asymmetries we 
are left with the 40 — 50% of the events. In both cases these ratios are good 
enough to allow the measurement of the observable of interest. Only in the case 
of variable invariant mass with fixed cut it's possible to distinguish the behavior 



of our model from that of the SSM (see Fig. 12 and compare with Fig. 11 1 
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Figure 4: Significance as a function of the corresponding cut associated to the 
four definitions of asymmetry for on-peak events in both the MiAUMSSM and 
the SSM (calculated with the Herwig package) 



4.4 Comparison with other models: LRM and SSM 

In this subsection we present a brief analysis of the results for the asymmetry 
obtained with two well-known models of extra U(l) extension of the SM: the 
Left-Right Model (LRM) PUCE] and the previously mentioned SSM [26]. We 
will see that the asymmetry in our anomalous model almost always leads to 
values which can be distinguished from those of the LRM and SSM. This implies 
that a possible future measure could discriminate among these models. 

The couplings among the fermions and the Z' for all these models could all 
be written in the form 

gz'J^ =aV^Vy-9i75)f/4 (19) 



where: 



2 cos t 



2 cos 0\Y 

g a MiAUMSSM 



SSM 

LRM (20) 



The charges gt and g f A are given explicitly in tablejij 6w is the Weinberg angle, 
defined by sin 6w = 0.231. In the case of the LRM we have chosen the so called 
symmetric version, for which aj,jj = 1.59 [TTJES. Using HERWIG, we have 
calculated the on-peak asymmetry for these two models. Obviously in the case 
of the MiAUMSSM we do not have a unique value for the asymmetry, because 
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Figure 5: Results for the asymmetries with Ql — for the best cuts. The left 
image represents Arfb, the right image represents Aq 
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Figure 7: Results for the asymmetries with Qq — for the best cuts. The left 
image represents Arfb, the right image represents Aq 
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Figure 9: Results for the asymmetries with Qh u = for the best cuts. The left 
image represents Arfb, the right image represents Aq 




Figure 10: Results for the asymmetries with Qh u = for the best cuts. The 
left image represents Ac, the right image represents A% 
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Figure 11: Ni/N to ti as a function of the corresponding cuts associated to the 
four definitions of asymmetry for on-peak events 
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Table 3: Couplings of the SM fermions to the Z's for the SSM, LRM and 
MiAUMSSM models. 



in the model the charges are not fixed. To show that it is possible to distinguish 
the MiAUMSSM from the other models we have to estimate the statistical error 
in this measurement. Using the formula (16), the statistical error is estimated 
for all the models we are considering. The exact values depend on the cross 
section which is model dependent. Now, if we fix Qh u = Qq = Ql = 0.5, the 
resulting values for the asymmetries associated to the three models are showed 
in figure [13] 

The data plotted in the figures show that it is always possible to discriminate 
the anomalous model from the non anomalous ones. 

Now we want to stress that the three charges of our model are free but 
the couplings of the fermions to the Z' in the anomalous MiAUMSSM have a 
peculiar functional form given in table [3] As a consequence it is not possible 
to match the couplings to the extra Z' of the MiAUMSSM with those of other 
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Figure 12: Ni/N toti as a function of the invariant mass, with the cuts kept fixed, 
for the four definitions of asymmetry 



models. But, since the four asymmetries have associated statistical errors we 
could have a range of values of our three charges where the couplings of the 
MiAUMSSM (and consequently the asymmetries) could be matched with those 
of the SSM and LRM models within the considered errors. In reality this does 
not happen as we can infer from Fig.s [14] and |15[ where we consider an error 
up to 25%, much bigger than the expected experimental error. Observing the 
amount of points in this figures, it is evident that the SSM is closer to our 
model than the LRM. This is the reason why troughout this paper we focus our 
analysis on the comparison with the SSM. Note that the colors in the images 
[14| and [15] are used only for guiding the eyes and are neither related to the 
asymmetry values nor express any functional relation among the charges. 



4.5 General case 

In this section we want to find the function which describes the asymmetry in 
terms of the three free charges of our model which can assume values between 
— 1 and 1. From the cross section of our process, that can be found in Appendix 
[A] we can see that the amplitude is proportional to the fourth power of the 
charges. So, the equations ([2]), (j3|), Q and Q imply that the asymmetry must 
be a rational function in which both the numerator and denominator are fourth 



14 



LRM| 



SSM | 



-0.4 -0.3 -0.2 -0.1 0.1 0.2 0.3 0.4 



MiAUMSSM | 



-0.1 0.1 0.2 0.3 



Figure 13: Left: Afb and Ao for the three models considered in the main text. 
Right:A c and A E . The charges of the MiAUMSSM are all fixed to 1/2. 





Figure 14: Values of the charges that give a MiAUMSSM asymmetry close to 
the SSM asymmetry within errors of 15, 20, 25 %. 
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Figure 15: Values of the charges that give a MiAUMSSM asymmetry close to 
the LRM asymmetry within errors of 20, 25 %. 



grade polynomials in the charges: 

T.l hk =^AQH u ) l (Q Q y{QL) k 

Ylij,k=o (Qh u y (Qq) j (QlY 1 j 

with i+j + k = n<4. 

The apparent symmetries of the contourplots obtained in subsection |4.2| 
imply that the terms of odd degree in the charges are suppressed. Then the only 
relevant terms which do not contain Qh u are Qq, Q\, Qq, Q l and QqQ\, while 
for the terms that do not contain Qq or Ql we can also have terms of the form 
QiQj or QiQj where i, j are the two free charges of each case. For example, a 
term proportional to QqQl is suppressed, while a term proportional to QqQh u 
is present. Fitting our data for the on peak asymmetries with the functional 



form (21), we find the coefficients a^'s of (21). Then considering only three of 
the four definitions we obtain a non-linear system with three equations and three 
variables (Qh u , Qq and Ql) which could be solved numerically. In this way the 
asymmetry is useful for fixing the values of the U(l)' charges. Moreover, once 
the values of the three charges are obtained by the previous system, the fourth 
definition of asymmetry can be used as a check for the validity of the model 
under exam. Infact its hypothetical experimental value must be recovered by 



using (21) with the values of the charges already found, within the considered 
error (we use the mean relative error (MRE) for each asymmetry definition). 
In the Appendix [B] we write a table with the coefficients of the four fits. As 
expected, we have found out that the odd degree polynomials have negligible 
coefficients, thus confirming the intuitions stemming from the analysis of the 
contour plots. The exactness of the fit is evaluated computing the R 2 [^] and 
the medium relative error for these results. The results are showed in table |4j 
attesting the accuracy of the procedure. In particular the R 2 value states (as we 
expected) that the errors in our fits are almost completely due to the numerical 
approximations in the calculation of the integrals and to the uncertainties of 
the PDFs. 



R is called coefficient of determination. A perfect fit has R = 1 
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Arfb 


Ao 


A c 


A E 


r 2 


0,999 


0,999 


0,999 


0,999 


AIRE 


0.008 


0.009 


0.019 


0.017 



Table 4: R 2 and Medium Relative Error for the polynomial fit of the asymmetry 
with respect to the three charges 

5 Conclusion 

We have numerically calculated the LHC asymmetry of the MiAUMSSM using 
four different definitions, namely forward-backward, one-side, central and edge 
asymmetries for the process pp e + e~ . An analogous study with a complete 
simulation of the ATLAS detector has been carried out for the SM asymmetry 
[27) . The performance of the detector will be unaltered for a measurement in 
the TeV range (hoping this to be the scale for the Z') meaning that our results 
should be robust even in the case of a complete analysis in which there is a 
full simulation of the detector. In this respect the good agreement, shown in 
Table 6.4 of [28], for the results coming from Pythia versus those coming from 
the real measurement testifies the good degree of precision of the Monte Carlo 
simulators for these kind of computations. In our case, to achieve a statistics 
comparable to that collected for the case studied in [37], at least 100 fb _1 will 
be needed and this will be achieved most likely in 2015 when the machine will 
work at */s = 14 - IhTeV after the planned shutdown in 2013/2014. 

To infer the optimal cuts to use we have maximized the significance related 
to each asymmetry. We have verified, as it is expected from [TU] , that these cuts 
are nearly independent from the free charges of our model. Then we have used 
these optimal cuts to investigate the asymmetry behaviour in function of pairs 
of free charges, keeping the third fixed to to have a graphical representation 
of the results. Furthermore we have found that the asymmetry is invariant 
under Qi —> —Qi- We have checked that even in presence of the cuts needed 
to evaluate the asymmetry, the number of left over events is such to lead to a 
meaningfull measurement. 

We have further shown that the MiAUMSSM is distinguishable from the 
SSM and LRM models, showing examples of different predictions for the asym- 
metries which differ at least for a good 20% of their value. Finally we have 
studied the four asymmetries as functions of the three charges and have fitted 
the results as a rational function of polynomials in the charges up to 4 th degree. 
The fit is found to be accurate, with a R 2 = 0.999 for all the definitions we have 
used. 
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A Cross section 

We have calculated the cross section in the CM for the process qq —> e + e~ for 
a general Drell-Yan interaction in which we can have the product of diagrams 



17 



q 



Figure 16: Feynmann diagrams for the Drell-Yan process in the MiAUMSSM 



in which the 7, Z and Z' can be exchanged . Thus we have 6 possible terms: 
77, "/Zq, jZ', ZqZq, ZqZ' and Z' Z' . The total amplitude is: 



\M\\ q ) 



1 1 

34 



2 dl 9b M ab 

a,b=~f,Z,Z' 



(22) 



where the fractions j and | come out from the averages over spin and color, 
go is the coupling associated to the Z' . We fix its value to 0.1. M ab is the 
amplitude of each process divided by the couplings: 



M. 



ab 
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(s - ml) (s - 
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T b m b ) 
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(r h m b ) 
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|2?7lg m 2 q (Cf~ a C eb — C^ a ^Y,b) ipq,a ^q,b ~ ^q,a^q,b) ~ m e (Pq ' Pq) x 
(pt,a C<^& — C<^b) (^a ^6 + ^q,a ^q.b) + (P« ' Pe) (Pq ' Pe) x 

(C& C e v Q + C e A a C e y b ) (C& C 9 y Q + C£ a Cl b ) + (p, • pe) (p 5 • p e ) X (23) 

(C<^a ^e.fc + C<Tb) (^«,a + ^-jS) _ (P« ' P e ) (P« ' P e ) X 



CA 
e,a ^ 



fiA r~iV \ 



(Pq ■ Pe) (Pq ' Pe) X 



C^a C<^b) (C^a C^fc + C<^a ~ TO <? ( Pe ' Ps ) X 

i°e,o °e,b "+" u e,a °e,&J l°q,a '-'g.fc — °g,o ^ q.b) ( 



In this expression N ab is a multiplicity factor that is equal to I if the exchanged 
vector bosons are identical and is equal to 1 if they are different. The C's are 
simply the vector and axial quantum numbers related to the vector bosons: for 
the 7 and the Z$ they are the usual SM quantum numbers that can be found in 
[2"5] , while the vector and axial couplings related to the Z' have been calculated 
in [53] and are showed in table [5j 

We remark that in this cross section there are only terms of fourth degree in 



the powers of the charges. However, from the equation (22 1 we know that there 
are different contributions to the total squared amplitude of our process. These 
terms are divided in three types: the term Z'Z', gives contribution of 4th degree 
in the anomalous charge; the terms ^Z' and ZqZ' give contributions of second 
degree in the anomalous charges; the terms 77, 7Z0 and ZqZq give contributions 
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f = u,c,t 


Qq + Q h J2 


-QhJ2 


f = d,s,b 


Qq-Q h J2 


QhJ2 


f = e,n,T 


Ql-QhJ2 


QhJ2 



Table 5: Vector and axial quantum numbers of the SM fermions with respect 
to the Z' 



of zero degree in the anomalous charges. Since we are studying the on-peak 
region, we expect the contribution from the Z'Z' channel to be dominant with 
respect to the others: this is evident from the table of the coefficients showed 
in the Appendix |Bj Observing the previous formula we can see that all the 
combinations of C's contain two C q and two C e because our elementary process 
involves two leptons and two quarks. Observing that Qq and Ql are related to 
C q and C e respectively, this implies that we cannot have terms of degree larger 
than 2 in Qq and Ql- This is verified by our fit, where the coefficients related 
to this terms are suppressed (see appendix [B]). The differential cross section can 
be found multiplying for the usual kinematic prefactor and summing this result 
over the contribution of the 6 possible initial quarks: 

s^L.'E^i'"''" < 24 » 

B Coefficients of the polynomial fit 

We have performed a numerical calculation of the asymmetries letting the three 
charges vary in the — 1 < Qi < 1 range. Then we have fitted the results with the 



rational function (21 ). We have found that only the even grade terms contribute 



to the results, so we neglect the odd grade terms. 



Another point to note is that the formula (21 1 implies that all the coefficients 
ciijk and bijk are defined up to a global multiplicative factor. To permit the 
comparison among the different types of asymmetries we have fixed a 40 o = 1 
(or a4oo = — 1 for the C asymmetry that assumes opposite sign with respect to 
the others). However, if such type of models will be discovered at the LHC, this 
value will be fixed differently to match the experimental results. The coefficients 
values for our choice are listed in tabic [U 
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Arfb 


Ao 


Ac 


A E 


aooo 


(-0.52 ±0.02) x 10~ B 


(-0.31 ±0.04) x 10~ 6 


(1.18 ±0.04) x lO" 6 


(0.86 ±0.18) x 10- e 


a 200 


(82 ± 3) x 10" 6 


(58 ±5) x 10" e 


(-17 ±5) x 10- e 


(140 ±21) x 10~ B 


O020 


(9.9 ± 1.5) x 10-" 


(9 ± 2) x 10" e 


(-27 ±3) x 10- e 


(12 ± 11) x 10" 6 


0(102 


(5.2 ± 1.4) x 10~ 6 


(2 ±2) x 10~ 6 


(11 ±2) x 10~ 6 


(61 ± 10) x 10"° 


dllO 


(5 ±3) x 10- B 


(4 ±5) x 10~ 6 


(19 ±6) x 10~ 6 


(113 ±24) x 10" 6 




(-18 ±4) x 10" B 


(-6 ±6) x 10" B 


(-148 ±7) x 10~ B 


(-269 ±27) x 10~ 6 


oon 


(-80 ±3) x 10- B 


(-65 ±5) x 10" 6 


(34 ±5) x 10~ 6 


(-177 ±22) x 10~ 6 


a 400 


1 (fixed) 


1 (fixed) 


-1 (fixed) 


1 (fixed) 


O040 


0.015638 ± 0.000004 


0.015634 ± 0.000006 


-0.015444 ±0.000007 


0.01552 ±0.00003 


0(104 


0.000969 ± 0.000002 


0.000967 ±0.000004 


-0.000984 ± 0.000004 


0.000964 ±0.000018 


O310 


0.88011 ± 0.00005 


0.87460 ± 0.00007 


-0.85003 ± 0.00008 


0.8573 ± 0.0003 


O220 


0.01525 ± 0.00004 


0.01544 ±0.00006 


-0.01565 ±0.00007 


0.0163 ± 0003 


Ol30 


-0.000729 ±0.000019 


-0.00061 ± 0.00003 


0.00077 ± 0.00003 


-0.00113 ±0.00013 


«301 


-1.99340 ±0.00005 


-1.99305 ± 0.00008 


1.99360 ± 0.00009 


-1.9930 ±0.0004 


a 2 n 


-1.75973 ±0.00010 


-1.74845 ± 0.00016 


1.69973 ± 0.00018 


-1.7141 ±0.0007 


0121 


-0.00388 ±0.00007 


-0.00403 ± 0.00011 


0.00429 ± 0.00012 


-0.0052 ± 0.0005 


O031 


0.00165 ± 0.00002 


0.00128 ±0.00003 


-0.00188 ± 0.00004 


0.00244 ±0.00015 


«202 


0.00456 ± 0.00010 


0.00438 ± 0.00016 


-0.00463 ± 0.00018 


0.0048 ± 0.0007 


0112 


-0.00049 ± 0.00009 


-0.00044 ± 0.00014 


0.00013 ± 0.00015 


-0.0011 ±0.0006 


O022 


0.00773 ± 0.00003 


0.07740 ± 0.00005 


-0.00778 ± 0.00006 


0.0068 ± 0.0002 


a 103 


-0.001244 ±0.000013 


-0.00120 ± 0.00002 


0.00136 ± 0.00002 


-0.00170 ±0.00009 


O013 


0.000181 ± 0.000015 


0.00023 ± 0.00002 


-0.00014 ±0.00003 


0.00047 ±0.00011 


&000 


(-1.19 ±0.06) x 10~ 6 


(-0.70 ±0.09) x 10~ 6 


(-3.19 ±0.12) x 10~ 6 


(2.2 ±0.4) x 10~ 6 


&200 


(121 ±7) x l(T e 


(181 ±12) x 10" e 


(-88 ± 16) x 10" e 


(-637 ±57) x 10- e 


&020 


(23 ±4) x 10~ 6 


(21 ±6) x 10~ 6 


(74 ± 7) x 10~ 6 


(29 ± 28) x 10" 6 


b()02 


(12 ±3) x 10~ 6 


(5 ±5) x 10- e 


(-29 ±6) x 10" B 


(154 ±25) x 10" e 


hio 


(-58 ±26) x 10~ 6 


(51 ±41) x 10" 6 


(-282 ±52) x 10~ 6 


(2392 ± 204) x 10~ 6 


hoi 


(-116 ± 13) x l(T e 


(-207 ±20) x 10" e 


(408 ± 27) x 10~ B 


(-392 ± 100) x 10" e 


ban 


(85 ± 39) x 10" b 


(-82 ±62) x 10~ 6 


(-104 ±79) x 10~ 6 


(-1120 ±306) x 10~ e 


&400 


1.90571 ± 0.00004 


1.90585 ±0.00006 


2.28541 ± 0.00008 


2.1465 ± 0.0003 


&040 


0.036021 ± 0.000010 


0.036021 ±0.000016 


0.04190 ± 0.00002 


0.03945 ± 0.00008 


&004 


0.002232 ± 0.000006 


0.002229 ± 0.000009 


0.002670 ± 0.000012 


0.00245 ± 0.00005 


&310 


1.40007 ± 0.00018 


1.3899 ±0.0003 


1.3347 ± 0.0004 


1.2650 ±0.0015 


&220 


3.8269 ± 0.0003 


3.8323 ± 0.0004 


4.5879 ± 0.0006 


4.338 ±0.002 


^130 


-0.00386 ±0.00018 


-0.0033 ± 0.0003 


-0.0035 ± 0.0004 


-0.0103 ±0.0014 


&301 


-3.79461 ± 0.00014 


-3.7943 ± 0.0002 


-4.5504 ±0.0003 


-4.2761 ±0.0011 


&211 


-2.7984 ±0.0004 


-2.7775 ± 0.0007 


-2.6661 ±0.0009 


-2.533 ±0.003 


&121 


-7.5882 ± 0.0007 


-7.5998 ±0.0011 


-9.1030 ±0.0014 


-8.592 ±0.005 


&031 


0.0081 ± 0.0002 


0.0061 ± 0.0003 


0.0098 ± 0.0005 


0.0168 ±0.0017 


&202 


3.8004 ± 0.0003 


3.8021 ± 0.0004 


4.5567 ± 0.0006 


4.291 ±0.002 


&112 


2.8003 ± 0.0005 


2.7871 ± 0.0009 


2.6705 ±0.0011 


2.524 ±0.04 


b()22 


7.6032 ± 0.0006 


7.6124 ±0.0009 


9.1181 ±0.0012 


8.599 ±0.004 


bl03 


-0.0023 ± 0.0002 


-0.0012 ± 0.0004 


-0.0026 ± 0.0005 


-0.0120 ±0.0018 


boi3 


-0.00009 ±0.00035 


0.0002 ± 0.006 


-0.0021 ± 0.0007 


0.017 ±0.003 



Table 6: Coefficients of the fits for the four definitions of asymmetry 
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Note that this table contains only the statistical error and not the systematic 
error due to the choice of the PDFs. 



References 



P. Langacker, "The Physics of Heavy Z' Gauge Bosons" , Rev. Mod. Phys. 
81 (2009) 1199 |[arXiv:0801. 1345 [hep^pE[ 



[1] 

[2] A. Lionetto and A. Racioppi, "Supersymmetry Breaking in a Minimal 

[3] 



Anomalous Extension of the MSSM", arXiv:1102.5040 hep-ph 



[4] 
[5] 
[6] 

[7] 
[8] 



P. Anastasopoulos, F. Fucito, A. Lionetto, G. Pradisi, A. Racioppi and 
Y. S. Stanev, "Minimal Anomalous U(l)-prime Extension of the MSSM", 
Phys. Rev. D 78 (2008) 085014 prXiv:0804.1156 [hep-th]].| 

F. Fucito, A. Lionetto, A. Mammareila and A. Racioppi, "Stueckelino 
dark matter in anomalous U(l)-prime models", Eur. Phys. J. C 69 (2010) 
455 |[arXiv:0811.1953 [hep-ph]]. | 



F. Fucito, A. Lionetto and A. Mammareila, "Dark Matter in Anoma- 
lous U{1)' Models with neutral mixing", Phys. Rev. D 84 (2011) 051702 



arXiv: 1105.4753 [hep-ph 



E. Komatsu et at [WMAP Collaboration], "Seven- Year Wilkinson Mi- 
crowave Anisotropy Probe (WMAP) Observations: Cosmological Inter- 
pretation", Astrophys. J. Suppl. 192 (2011) 18 |[arXiv: 1001 .4538 [astn>| 
|ph.CO]].| 

A. Lionetto and A. Racioppi, "Gaugino radiative decay in an anomalous 
U(l)- prime model", Nucl. Phys. B 831 (2010) 329 |[arXiv:0905.4607 [hep^| 

F. Fucito, A. Lionetto, A. Racioppi and D. R. Pacifici, "A Phenomenolog- 
ical study on the wino radiative decay in anomalous U(l)' models", Phys. 



Rev. D 82 (2010) 115004 [arXiv:1007.5443 [hep-ph 



[9] P. Langacker, R. W. Robinett, J. L. Rosner, "New Heavy Gauge Bosons 



in p p and p anti-p Collisions" , |Phys. Rev. D30 (1984) 1470. 



[10' 



Z. q. Zhou, B. Xiao, Y. k. Wang and S. h. Zhu, "Discriminating Differ- 
ent Z's via Asymmetries a t the LHC" , Phys. Rev. D 83 (2011) 094022 
|[arXiv: 1102. 1044 [hep-ph"]J] 



[11] F. Petriello and S. Quackenbush, "Measuring Z' couplings at the CERN 



LHC", Phys. Rev. D 77 (2008) 115004 [arXiv:0801.4389 [hep-ph 



[12] M. Cvetic and S. Godfrey, "Discovery and identification of extra gauge 



bosons" , |arXiv:hep -ph/9504216. 



[13] 



M. Dittmar, A. S. Nicollerat and A. Djouadi, "Z-prime studies at the LHC: 
An Update", Phys. Lett. B 583 (2004) 111 |[arXiv:hep-ph/0307020]Tl 



21 



[14] S. Godfrey, T. A. W. Martin, "Identification of Extra Neutral Gauge 
Bosons at the LHC Using b- and t-Q uarks" , Phys. Rev. Lett. 101, 151803 
(2008). [|arXiv:0807. 1080 [hep-ph]]T] 

[15] Y. -k. Wang, B. Xiao, S. -h. Zhu, "One-side forward-backward asymmetry 



at the LHC", Phys. Rev. D83 (2011) 015002. [arXiv:1011.1428 [hep-pli 



[16] Y. -k. Wang, B. Xiao, S. -h. Zhu, "One-side Forward-backward Asym- 
metry in Top Quark Pair Production at CERN Large Hadron Collider", 
Phys. Rev. D82 (2010) 094011. [jarXiv: 1008.2685 [hep-phj].| 

[17] P. Ferrario and G. Rodrigo, "Charge asymmetries of top quarks: A Win- 
dow to new physics at hadron colliders", J. Phys. Conf. Ser. 171 (2009) 
012091 1" [arXiv:0907.0096 [hep-ph]]^ 



[18] J. H. Kuhn, G. Rodrigo, "Charge asymmetry in hadroproduction of heavy 



quarks", Phys. Rev. Lett. 81 (1998) 49-52. [arXiv: hep-ph/9802268 



[19] J. H. Kuhn, G. Rodrigo, "Charge asymmetry of heavy quarks at hadron 
colliders", Phys. Rev. D59 (1999) 054017. [[hep-ph/9807420]T1 

[20] O. Antunano, J. H. Kuhn, G. Rodrigo, "Top quarks, axigluons and 
charge asym metries at hadron colliders", Phys. Rev. D77 (2008) 014003. 
|[arXiv:0709.1652 [hep-phjj.| 

[21] P. Ferrario, G. Rodrigo, "Massive color-octet bosons and the charge asym- 
metries of top quarks at had ron colliders" , Phys. Rev. D78 (2008) 094018. 
|[arXiv:0809.3354 [hep-ph]].| 

[22] B. Xiao, Y. K. Wang, Z. Q. Zhou and S. h. Zhu, "Edge Charge Asymmetry 
in Top Pair Production at the LHC", Phys. Rev. D 83 (2011) 057503 



[arX iv:1101.2507 [hep-ph]] . 



[23] A. Mammarella, Ph.D. Thesis, "Anomalous U(l), Dark Matter and Asym- 
metry" , i arXiv: 1201 .4340 [hep-phfl 

[24] G. Corcella, I. G. Knowles, G. Marchesini, S. Moretti, K. Odagiri, 
P. Richardson, M. H. Seymour and B. R. Webber, "HERWIG 6: An Event 
generator for hadron emission reactions with interfering gluons (including 



supersymmetric processes)", JHEP 0101 (2001) 010 [hep-ph/0011363] 



[25] G. Corcella, I. G. Knowles, G. Marchesini, S. Moretti, K. Odagiri, 
P. Richardson, M. H. Seymour and B. R. Webber, "HERWIG 6.5 release 



note", hep-ph/0210213. 



[26] G. Altarelli, B. Mele and M. Ruiz-Altaba, "Searching for new heavy vector 
bosons in p anti-p colliders" , |Z. Phys. C 45 (1989) 109| [Erratum-ibid. C 
47 (1990) 676]. 



[27] G. Aad et al. [ATLAS Collaboration], arXiv: 1203.3100 [hep-ex 



[28] Cattani, G., PhD Thesis, "Forward-backward asymmetry measurement in 
pp — » Z/j* — > /J + /i~ events at ^/s = 7TeV with the ATLAS experiment 



at the LHC," CERN-THESIS-2011-247. 



22 



[29] F. Halzen, A.D. Martin, "Quark and Leptons: an Introductory Course to 
Modern Particle Physics" (1984) 



23 



